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0.  Abstract.  R  and  S  charts  for  monitoring  the  variability  of  a 

2 

process  are  studied  when  the  control  variance,  a^,  is  unknown,  when 

observations  are  non-Gaussian,  and  when  the  process  drifts.  Geometric 

2 

bounds  on  run-length  distributions  for  the  R  and  S  charts  are  given  for 

2 

all  underlying  distributions  when  is  estimated  in  a  base  period. 

Similar  bounds  apply  in  the  case  of  a  drifting  process,  the  notion  of 

which  is  developed  rigorously.  It  is  shown  that  normal-theory  properties 

of  these  charts  hold  exactly  for  every  spherical  process  when  a  suitable 
2 

estimate  for  is  used.  Run-length  distributions  otherwise  are  shown 
to  be  ordered  stochastically  given  a  peakedness  ordering  on  the  under¬ 
lying  spherical  processes. 
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1.  Introduction.  Standard  procedures  for  monitoring  the  variability 

2 

of  a  production  process  are  the  R  and  S  charts.  Let  {Xi  «  (X^,  X  2,  .... 

XiQ);  i  -  1,  2,  ...}  be  the  outcomes  and  {[X^^  —  Xi(2>  —  *’*  —  Xi(n)^’ 

i  *  ^  2,  ...}  the  ordered  outcomes  in  a  succession  of  samples  of  n 

2  2 

observations  each;  let  o  be  the  process  variance  and  Oq  its  value  when 

the  process  is  in  control;  and  let  ■  ^Xi(n)  ”  ^i(l)^  an<*  ^i  “  ^j«l^Xij  “ 

2 

3^)  /(n-1)  respectively  be  the  typical  sample  range  and  sample  variance. 

2 

The  R  and  S  charts  are  graphs  of  the  standardized  ratios  {X^A^q*  i  ■  1,  2, 

2  2 

...}  and  (S^/cig;  i  ■  1,  2,  ...}  plotted  against  time  on  the  horizontal 

scale.  The  monitored  process  is  taken  to  be  in  control  as  long  as  the 

2  2 

ratios  are  within  their  control  limits,  i.e.,  R ./an  <  c  or  S./a_  <  d  ,  as 

lu  —  ct  i  0  —  a 

appropriate.  Otherwise  the  chart  signals  that  the  process  is  not  in  con¬ 
trol  and  that  corrective  action  is  needed.  The  object  is  to  detect  early 
a  shift  in  the  variance,  yet  to  signal  infrequently  when  the  process  is 
in  control.  The  sample  number  N  at  which  a  signal  occurs  is  called  the  run 

length;  its  distribution  determines  the  operating  characteristics  of  the  R 
2 

and  S  charts.  The  mean  of  the  run- length  distribution,  if  defined,  is 
called  the  average  run  length. 

2 

Usually  it  is  assumed  that  Oq  is  known;  that  independent  random  samples 

are  generated  on  successive  occasions  by  a  Gaussian  process;  and,  for 

2 

evaluating  run-length  distributions,  that  a  remains  constant  throughout 

'l 

the  monitoring  period.  Then  the  run-length  distributions  for  the  R  and 

charts  are  geometric,  their  parameters  being  the  probability  of  exceeding 

2 

the  control  limit  on  any  particular  occasion.  Although  S  is  more  efficient 
than  R,  the  latter  claims  wide  usage  in  statistical  quality  control  owing  to 
its  greater  simplicity  and  the  fact  that  loss  in  efficiency  may  be  compensated 
by  larger  samples  in  nondestructive  testing. 

2 

In  routine  use  one  or  more  of  the  usual  assumptions  may  fail.  When 


3 


is  unknown  a  common  practice  is  to  estimate  it  during  a  base-line  period 

2 

when  the  process  is  in  control,  and  to  maintain  the  R  and  S  charts  using 

2 

this  estimate  in  place  of  o^.  The  run-length  distributions,  no  longer 

geometric,  are  complicated  by  stochastic  dependence  among  the  resulting 

ratios.  Properties  of  the  actual  run-length  distributions  are  examined 

2 

in  Section  2  using  a  variety  of  estimates  for  a  without  stipulating  the 

0 

2 

underlying  distribution.  Even  when  Oq  is  known  the  run-length  distribu- 

2 

tions  are  not  geometric  unless  ?  remains  constant.  The  notion  of  a 

drifting  process  is  developed  precisely  in  Section  2  and  stochastic  bounds 

for  its  run-length  distributions  are  given.  When  Gaussian  assumptions  fail 

little  is  known  beyond  the  usual  caution  that  normal-theory  procedures  for 

variances  are  not  robust.  In  Section  3  the  robustness  of  normal-theory 
2 

R  and  S  charts  is  studied  for  a  large  class  of  underlying  processes. 

Section  4  gives  some  numerical  comparisons  of  practical  interest  stemming 

2 

from  earlier  sections,  with  special  reference  to  normal-theory  R  and  S  charts. 

2.  Bounds  on  Run-length  Distributions.  Given  a  sample  *om  (xoi-* 

X^)  of  m  observations  taken  while  the  process  is  in  control,  let  be  an 

estimate  for  a^.  In  particular,  may  be  the  sample  range  R^,  the  sample 

standard  deviation  S^,  the  inter-quartile  range,  or  some  other  estimate 

which  need  not  be  stipulated  here.  Denote  by  N  and  Nc  the  run  lengths  of 
2 

the  R  and  S  charts,  and  by  F  (t)  and  F_(t)  their  respective  cumulative 

J\  J 

distribution  functions  (cdf's).  When  the  process  variance  remains  constant 
throughout  the  monitoring  period  these  distributions  are  determined  by  the 
Studentlzed  ratios  {R^/V^;  i  ■  1,  2,  ...}  and  i  *  1,  2,  . . . }  as 

follows,  where  P(Ejv)  and  P(E)  are  conditional  and  unconditional  probabili¬ 
ties  of  the  event  E.  It  is  clear  for  each  positive  integer  t  that 
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P(N  >  t)  -  P(R./V  <  c  ,  ....  R  /V  <  c  ) 

R  e,10  t  0  -  a  (2>1) 

-  /"  [P(R  <  vc  |v)]tdG(v) 

0  01 

P(N.  >  t)  -  P(S?/V2  <  d  ,  ....  S 2/V2  <  d  ) 

S  10—  a  t  0  —  a 

(2.2) 

-  /I  [PCS2  <  V2d  | v) ]tdG(v) 

0  —  a 

2  2 
where  R  and  S  have  the  same  distributions  as  R^  and  S^,  and  G(*)  is  the 

cdf  of  VQ.  Moreover,  in  support  of  later  developments  using  finite-dimen¬ 
sional  distribution  theory,  we  note  that  P(N  <  =°)  ■  1  *  P(N  <  <*>) ,  i.e., 

K.  b 

both  charts  eventually  signal  with  unit  probability.  For  R  charts  this 
follows  on  writing 

P(Nn  »  t)  «  £  [P(R  <  vc  |v)]t-1{P(R  >  vclv)]  dG(v);  (2.3) 

iv  u  ct  at 

on  noting  that  H(t;v)  »  [P(R  <  vc  |v)]C  ^[P(R  >  vc  | v) ]  <  1,  that  /  dG(v)  ■  1, 
and  that  H(t;v)  ->  0  as  t  -*■  00  for  each  v  >_  0;  and  on  applying  the  dominated 
convergence  theorem  to  infer  that 

lim  P (N  -  t)  »  fl  lim  H(t;v)dG(v)  -  0.  (2.4) 

t-^oo  ^  ^  t  00 

Similar  arguments  establish  that  P(N  <  «)  ■  1. 

S 

Under  Gaussian  assumptions  a  survey  of  the  relevant  distribution  theory 
follows.  With  V-  ■  Sn  the  probability  P(N  >  t)  is  determined  by  the  joint 
distribution  of  Studentized  ranges  whose  marginals  are  of  the  type  considered 
and  tabulated  by  Pachares  (1959),  Harter  (1960),  and  Owen  (1962),  although 
the  joint  distribution  appears  not  to  have  been  studied.  Similarly  P(Ng  >  t) 
derives  from  the  known  Studentized  largest  chi-squared  distribution  considered 
by  Armitage  and  Kris hnaiah  (1964),  David  (1956),  Finney  (1941),  Ghosh  (1955), 


Gupta  (1963),  Nair  (1948),  and  Ramachandran  (1956).  Under  Gaussian 

assumptions  with  Vn  ■  R  ,  P(N  >  t)  is  determined  by  the  joint  distri- 

bution  of  (R^/Rg,  ...,  R^R^l ,  apparently  not  considered  before,  each 

marginal  of  which  is  of  the  type  studied  by  Link  (1950) .  To  evaluate 

(2.1)  and  (2.2)  in  particular  cases,  recourse  may  be  taken  to  numerical 

methods;  unfortunately,  even  the  extensive  tables  of  Armitage  and  Krishnaiah 

(1964)  are  tabulations  of  upper  percentiles  limited  to  the  range  1  <_  t  <_  12 

and  thus  are  not  suited  to  evaluating  P(N<,  >  t)  . 

Rather  than  treat  normal- theory  distributions  further,  we  turn  instead 

to  general  structural  properties  of  F  (t)  and  F  (t)  which  depend  neither 

R  b 

on  Gaussian  assumptions  nor  the  particular  choice  for  V^.  To  this  end  let 

G(t;a)  be  the  cdf  of  the  geometric  distribution  having  the  parameter  a.  Our 

first  result  provides  bounds  on  the  actual  run-length  distributions  F  (t) 

R 

and  Fg (t)  in  terms  of  G(t;a). 

THEOREM  2.1.  Let  {X^,  X^,  . « . ,  X  }  be  mutually  independent  sample  outcomes 

such  that  (X,,  ...,  X  }  are  distributed  identically.  Choose  c  and  d  such 
~1  a  a 

2  2 

that  P(R. /V„  <c)»l-a=  P(S,/V„  <  d  ),  and  let  N  have  the  distribution 
i  0  —  a  i  0  —  a 

G(t;a).  Then  N  and  N  are  stochastically  larger  then  N  whatever  be  the 
R  S 

underlying  distribution  and  choice  for  V  ,  i.e.,  G(t;a)  >  F  (t)  and  G(t;a) 

U  K 

>  F  (t)  for  every  t  >  0. 

b 

2  2 

Proof .  Consider  sequences  i  a  1,  2,  ...}  and  (S^/V^;  i  *  1,  2,  ... 

2  2 

of  independent  ratios  having  the  same  distributions  as  an<*  $^/Vq, 

*  * 

respectively,  and  let  and  be  run  lengths  associated  with  these. 

*  * 

Because  and  clearly  are  geometric,  it  suffices  to  show  that  >  t) 

*  * 

>  P(Nr  >  t)  and  P(Ng  >  t)  >  P(Ng  >  t)  for  each  fixed  t.  Consider  some 
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function  Xg)  >_  0  and  suppose  X  has  the  distribution  of  X^.  Under 

the  hypothesis  of  the  theorem  it  is  known  (cf .  Khatri  (1967),  for  example) 
that 

PWXj.Xg)  <  c,  ....  4>(Xc,X0)  <  c)  >_  [P (4> (X,Xg)  <_  c]\  (2.5) 

2  2 

On  choosing  4>  successively  as  $(X^,Xg)  =*  R^/Vg  ant*  ‘K^.Xq)  =  S^/Vg  ancl 

applying  (2.5)  to  (2.1)  and  (2.2),  we  find  that  P(NR  >  t)  >_  [P(R/Vq  <_  c^)  ] C 

*  * 

*  P(N  >  t)  and  P(N_  >  t)  >  P(N„  >  t) ,  as  required. 

R  b  —  b 

Several  consequences  of  Theorem  2.1  are  immediate.  One  is  that  average 

2  2 

run  lengths  for  the  R  and  S  charts,  when  cjg  is  estimated,  are  greater  than 

those  of  the  approximating  geometric  distributions.  Of  considerable  practical 

interest  is  the  error  incurred  on  approximating  FR(t)  and  Fg(t)  by  G(t;a); 

this  topic  is  studied  numerically  in  Section  4  for  the  Gaussian  case. 

To  make  further  progress  we  suppose  the  underlying  distributions  belong 

2  2 

to  a  scale-parameter  family.  Let  {cjg,  a^,  ...}  be  a  sequence  of  scale 

2,  2 

parameters  corresponding  to  successive  samples  and,  with  y^  =  cr^/ag,  let 

y  *  (Yi,y2,  ...).  Recall  that  expressions  (2.1),  (2.2)  and  Theorem  2.1  apply 

when  {X^,  ...,  Xfc}  are  distributed  identically;  under  a  scale-parameter 

2  2 

family  this  corresponds  either  to  the  case  =  a^,  when  the  process  is  in 

2  2  2  2 

control,  or  to  the  case  •  a  for  all  i  and  some  o  >  To  treat  the 

general  case,  here  called  a  drifting  process ,  let  T  consist  of  all  sequences 

Y  =  (Yj>  Y2»  •••)  for  which  y^  ^  1  for  all  i.  Two  such  sequences  are  said 

to  be  ordered  as  y  y'  whenever  y  <_  y^  for  all  i.  Under  this  incomplete 

ordering  an  element  y  e  T„  is  said  to  be  minimal  for  r„  c  r.  and  y..  e  T_  to 

-m  u  -  u  -MO 

be  maximal  for  T.,  if,  for  every  y  e  r  ,  y_-<  y y  .  In  particular,  y  = 

w  —  u  -U*  ^  _  in 

(1,  1,  ...)  is  minimal  for  T. 

As  the  run-length  distributions  for  drifting  processes  generally  depend 


L. 
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on  y ,  we  write  the  typical  cdf  F(t)  as  F(t;y)  and  consider  the  family 

(F(t;y);  Y  £  F}.  Such  a  family  is  said  to  be  stochastically  decreasing 

in  Y  if,  for  Y  and  Y'  in  T,  the  ordering  Y  -<  Y'  implies  the  stochastic 

ordering  F(t;  Y")  >  F(t;Y)  for  every  t  >  0.  That  the  families  (F  (t;Y); 

Y  £  F}  and  {Fg(t;Y);  Y  £  F}  of  run-length  distributions  are  stochastically 

2 

decreasing  in  Y,  and  thus  that  the  R  and  S  charts  tend  to  signal  more 
frequently  as  the  process  drifts  further  from  control,  is  established  in 
the  following  theorem  for  scale-parameter  families  of  distributions. 


THEOREM  2.2.  Let  {F„(t;y);Y  e  F}  and  {F„(t;y);Y  e  F}  be  families  of  run- 

R  ^  ^  j  «  w 

2 

length  distributions  for  the  R  and  S  charts  when  (X^,  ,  ...}  are  in¬ 

dependent  samples  from  a  scale-parameter  family.  Then  these  run-length 
distributions  are  stochastically  decreasing  in  y- 

Proof.  Consider  F„(t;Y),  the  arguments  for  F„(t;Y)  being  identical. 

R  O 

Corresponding  to  (2.1)  let  F^(u^,  ...,  ut>  be  the  joint  cdf  of  Cr^/V^,..., 

Rt/Vg}  in  the  general  case;  let  F^(u^,  ,  ut>  be  the  case  that  Y^  =  1 

for  all  i;  and  write  P  (NR  >  t)  =  F^(ca,  ...,  ca) .  By  independence  and 
the  assumption  of  a  scale-parameter  family,  we  infer  that 


Vnr  >  c) 


(2.6) 


The  proof  is  completed  on  noting  that  P  (N  >  t)  is  a  decreasing  function 

y  R 


of  (y^  Y2> 


.)  and  thus  F  (t;Y)  is  an  increasing  function,  as  asserted. 
R  — 


Run- length  distributions,  depending  on  an  excess  of  parameters,  are 

2 

complicated  in  the  case  of  drifting  processes  even  when  is  known.  Under 

further  conditions  on  Y  Theorem  2.2  supports  the  construction  of  envelopes 

of  curves  containing  F  (t;Y)  and  F  (t;Y)  in  terms  of  distributions  of  types 

R  5 

(2.1)  and  (2.2).  Suppose  Y  =  (Y ^ ,  Y^»  •••)  is  a  bounded  sequence  having 
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Y  =  inf (y . ,  Y „t  •••)  and  y  -  sup(Y. ,  Y0,  •  ••),  and  let  y(a)  «  y  (1,  1,  ...) 
mii  Miz  -  m 

and  y(M)  =  ym(1,  !•»  •••).  Clearly  FR(t;y(m))  and  FR(t ;y (M) )  are  distributions 

of  type  (2.1),  and  Fg(t;  y(o))  and  Fg(t;y(M))  are  of  type  (2.2),  these  being 
2 

geometric  when  an  is  known.  Bounds  on  F  (t;y)  and  F  (t;y)  for  drifting 
U  R  -  S  ~ 

processes  are  given  explicitly  in  the  following  theorem  for  any  scale- 
parameter  family  of  underlying  distributions. 


THEOREM  2.3.  Let  F„(t;y)  and  F„(t;Y)  be  run-length  distributions  for  the 
R  S'* 

2 

R  and  S  charts  of  a  drifting  process  having  the  bounded  parameter  sequence 

Y  =  (y^.  Y2»  •■•)•  If  (Xq,  Xl,  ...}  are  independent  samples  from  a  scale- 

parameter  family,  then  both  F  (t;y)  and  F  (t;y)  satisfy  inequalities  of  the 

R  b  ~ 

type 

(i)  F ( t ; y (M) )  >  F(t;Y)  >  F(t;y(m)) 

2 

for  every  t  >  0.  In  particular,  when  a ^  is  known  it  is  true  that 

(ii)  G ( t ; a (M) )  >  F(t;y)  >  G(t;o(m)) 

where  a(m)  and  ct(M)  respectively  are  the  probabilities  of  exceeding  the 
control  limits  on  any  particular  occasion  when  y^  =  Ym  and  when  y^  =  Y^  for 
all  i. 

Proof .  The  theorem  follows  directly  from  the  monotonicity  properties  of 

F  (t;y)  and  F  (t;y)  established  in  Theorem  2.2,  together  with  the  ordering 
R  b  ~ 

Y(m)  U.  y«<y(M)  under  the  condition  that  y  be  bounded. 


Several  conclusions  follow  immediately.  Because  Yq 


(1,  1,  •••)  is 


minimal  in  F,  it  follows  for  scale-parameter  families  that  P  (NR  <  ®)  =  1  = 

2 

P  (N„  <  ®)  and  thus  the  R  and  S  charts  eventually  signal  with  unit  probability 
Y  « 

whatever  be  the  parameter  sequence  in  T.  Moreover,  the  lower  and  upper  bounds 
of  Theorem  2.3  become  tighter  as  the  process  drifts  between  narrowing  limits, 
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i.e.,  as  y  -  y  decreases.  The  extent  of  this  tightening  depends  on 
Mm  * 

the  underlying  distributions  and  is  studied  numerically  in  Section  4 

2 

with  reference  to  normal-theory  R  and  S  charts.  We  finally  note  that, 
if  Tq  consists  of  all  sequences  y  for  which  y(m)  is  minimal  and  y(M)  is 
maximal,  then  the  given  bounds  apply  uniformly  to  all  members  of  the 
family  {F(t;y);  y  e  rQ}. 

2  2 

3.  Robust  Aspects  of  R  and  S  Charts.  In  practice  R  and  S  charts 

typically  are  based  on  Gaussian  assumptions  together  with  the  independence 

of  {Xq,  X^,  ...}.  Often  neither  assumption  is  tenable  in  sampling  from  a 

continuing  production  process.  Here  the  weaker  assumption  of  spherical 

symmetry  is  made,  allowing  at  once  dependencies  among  observations  and  a 

large  class  of  alternative  distributions.  We  suppose  first  that  successive 

observations  are  generated  by  a  spherical  process  and  thus  are  dependent 

in  all  except  the  Gaussian  case.  An  alternative  model,  for  use  when 

successive  samples  are  known  to  be  independent,  replaces  the  normality  of 

each  sample  of  size  n  by  n-dimensional  spherical  symmetry. 

A  distribution  L(X)  on  the  N-diraensional  Euclidean  space  R1*  is  said 

to  be  spherical  if,  for  each  (N  x  N)  orthogonal  matrix  Q,  L(QX)  =  L(X) . 

Each  such  X  admits  the  representation  X  =»  RZ  in  polar  form  such  that  L(X) 

is  determined  completely  by  the  radial  distribution  L(R)  on  (0,®) ,  and  Z 

N 

is  distributed  uniformly  on  the  surface  of  the  unit  sphere  in  R  independently 

of  R.  Let  SN(0,IN)  be  the  class  of  N-dimensional  spherical  laws  symmetric 

about  9  e  RN,  i.e.,  if  L( Y)  e  S^CS.I^),  then  L(Y  -  8)  is  spherical;  let  M 

m  N  k 

be  a  linear  subspace  of  k  ;  and  call  a  function  <p:  R  -*■  R  as  location-invariant 
with  respect  to  M  if,  for  every  5  e  M  and  X  e  R^,  MX  +  O  =  MX).  The 
following  result  is  basic  to  our  subsequent  developments. 
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LEMMA  3.1.  Let  L( Y)  be  spherically  symmetric  about  0  e  M,  and  let 
N  k 

ip:  R  -*■  R  be  any  location- invariant  function  homogeneous  of  degree 

k 

zero.  Then  the  distribution  L(p(Y))  on  R  is  invariant  for  all  distri¬ 
butions  L(Y)  in  the  class  {S„(9,I„);  9  e  M}. 

Proof .  By  the  location- invariance  of  d(')  and  the  fact  that  M  is  a 
subspace  of  R^,  we  may  assume  that  9  »  0.  Then  J  admits  the  polar 
representation  Y  *  RZ.  We  have  L(<j>00)  *  L($(RZ))  and,  by  homogeneity, 
L(<J>(RZ))  =  L(p(Z)).  The  fact  that  Z  is  uniform  on  the  surface  of  the  unit 
sphere,  independently  of  the  particular  L(Y)  e  S^CQ.I^),  assures  invariance 
of  the  distribution  L(<( >(Y))  as  asserted. 

Let  N  be  countable  and  let  (U(t);  t  e  M)  be  a  stochastic  process  such 
that,  for  each  N  and  {t^,  ...»  t^}  e  W,  the  joint  distribution  of  [U(t^), 


...,  U C t N) ]  is  spherical;  then  (U(t);  t  e  SI}  is  called  a  spherical  process . 
It  is  known  (cf.  Hartman  and  Wintner  (1940),  for  example)  that  a  countable 
process  is  spherical  if  and  only  if  it  is  a  scale  mixture  of  spherical 
Gaussian  processes,  i.e.,  if  for  each  N  and  {t^,  •••>  t^}  e  M,  the  joint 
probability  density  function  (pdf)  of  [U(t^),  ....  U(tN)]  is  given  by 


f(ur  ....  uN)  -  (2u)'N/2  /”T_Ne'-'-/2t2dG(T) 


(3.1) 


(u,  ,  ...,  u  )  and  G(-)  is  some  cdf  on  (0,°°).  Various  choices 
I  N 


where  u 

for  G(*)  yield  as  special  cases  the  spherical  Gaussian  process,  a  spherical 

Cauchy  process  whose  marginals  are  spherical  Cauchy  laws,  and  other  spherical 

stable  processes  as  further  examples. 

Against  this  background  we  consider  the  performance  of  normal-theory  R 
2 

and  S  charts  when,  in  fact,  the  distribution  of  {Xq,  Xj  ,  ...}  exhibits 
spherical  symmetry.  Partition  W  into  the  disjoint  subsets  (Nq,  . ..} 
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corresponding  to  successive  sampling  periods;  let  y(t)  and  o(t)  be 
functions  on  W  such  that  o(t)  >  0;  and  model  the  process  {X(t);  t  e  W} 
as  X(t)  *  o(t)U(t)  +  p(t)  such  that  (U(t);  t  e  W}  is  a  spherical  pro¬ 
cess.  We  assume  that  p(t)  and  cr(t)  are  step  functions,  constant  within 
samples  with  jumps  occurring  between  samples,  as  it  might  happen  when 
periods  of  intense  sampling  activity  are  followed  by  long  interim  periods 
and  the  functions  are  slowly  varying.  Accordingly  let  *  y(t)  and  o  * 

o(t)  for  t  e  i  *  0,  1,  2,  ...»  and  let  y  »  (y^,  •••)  such  that 

2  2  2 
y^  =  o^/Oq  >.  1.  Suppose  for  the  present  that  is  not  known  and  that 

2 

some  estimate  V  from  a  base-line  period  is  used,  and  denote  by  G  (t;y) 
u  R 

2 

and  Gg(t;y)  the  actual  run-length  distributions  of  the  R  and  S  charts 
under  Gaussian  assumptions.  Under  these  conditions  normal-theory  pro¬ 
cedures  are  more  than  robust;  they  remain  exact  for  all  spherical  processes 
as  shown  in  the  theorem  following. 

THEOREM  3.1.  Suppose  {XQ,  X^,  ...}  are  generated  from  a  spherical  process 

{U ( t) ;  t  e  N}  as  {X±(t)  «  a±U(t)  +  ji  ;  t  e  for  i  =  0,  1,  2 .  If 

Oq  is  estimated  using  a  location-invariant  function  of  X^,  then  the  run- 

2 

length  distributions  of  the  R  and  S  charts  are  identical  to  their  normal- 
theory  forms  G  (t;y)  and  G  (t;y),  respectively,  whatever  be  the  underlying 
spherical  process  (U(t);  t  e  N). 

Proof.  For  fixed  t  let  N  *  m  +  tn  be  the  number  of  observations  taken  up 

to  the  point  that  a  signal  occurs,  and  let  Y'  =  (Y^,  Y^,  ...»  Y^)  such  that 

-1  N 

Y^  *  X^.  Then  L(Y)  e  SN<8,IN)  for  0  belonging  to  a  subspace  M  of  R  , 

and  observe  that  ij>(Y)  *  C ,  ....  y  *Rt/V  )  not  only  is  location- 

invariant  with  respect  to  M,  but  is  homogeneous  of  degree  zero.  From  Lemma 

3.1  it  follows  that  the  joint  distribution  of  ...»  yt^Rt/Vg),  and 

thus  of  (R^/Vq,  ....  R^/Vq) ,  is  invariant  for  all  L(Y)  e  S^(9,I^).  The  joint 
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distribution  of  (R^/V^,  ....  Rt/Vg)  thus  is  independent  of  the  particular 

spherical  process  generating  the  observations  and,  as  the  spherical 

Gaussian  process  belongs  to  the  class,  the  run- length  distribution  P  (N  <  t) 

y  R 

is  identical  to  its  normal-theory  form,  namely,  G„(t;y)-  Similar  arguments 

2 

establish  invariance  properties  of  S  charts. 

It  may  be  noted  that  normal-theory  properties  of  the  run-length 

distributions  G  (t;y)  and  G  (t;y)  were  established  in  the  preceding  section 

R  S  **' 

under  the  assumption  that  (x^,  X^,  . ..}  are  mutually  independent.  Theorem  3.1 
assures  that  these  properties  continue  to  hold  for  all  spherical  processes 
despite  dependencies  among  successive  samples,  as  long  as  from  a  base¬ 
line  period  is  used  to  estimate  cr^.  In  contrast  to  these  results,  the 

2 

run-length  distributions  of  the  R  and  S  charts  do  depend  on  the  underlying 
spherical  process  when  cr^  is  known  and  the  base-line  period  thus  remains 
unobserved.  Despite  dependencies  among  successive  samples,  stochastic 
bounds  for  run-length  distributions  analogous  to  those  of  Theorem  2.1  are 
given  as  follows  for  spherical  processes. 

THEOREM  3.2.  Suppose  {X^,  X^,  ...}  are  generated  from  a  spherical  process 
2  2 

for  which  Oq  is  known  and  a  remains  constant,  and  let  NR  and  be  run 

2 

lengths  for  the  R  and  S  charts  having  the  respective  cdf's  F  (t)  and  F  (t) . 

K  b 

2  2  2 

Choose  ■  P(R/a_  >  c  )  and  ae  ■  P(S  /an  >  d  ),  where  R  and  S  have  the 

K  U  Cl  b  0  Cl 

2  it  it 

same  distributions  as  R^  and  for  i  *  1,  2,  ...  .  Let  NR  and  Ng  have  the 

geometric  distributions  G(t;  aB)  and  G(t;  ae).  Then  N  and  N  respectively 

K  b  Kb 

*  * 

are  stochastically  larger  than  N  and  N  ,  i.e., 

K  S 

G(t;aR)  >.FR(t),  G(t;  ag)  >_  Fg(t) 


for  all  t  >  0. 
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Proof.  We  supply  details  for  R  charts;  the  proof  for  S  charts  is 
identical  using  appropriate  notation.  On  representing  a  spherical  process 
as  a  scale  mixture  of  Gaussian  processes,  we  infer  that  the  actual  run- 
length  distribution  can  be  expressed  as  the  scale  mixture 

P(Nr  >  t)  -  /*  PCRj/tOq  <_  c^ . Rt^T<50  ±cctlT)dG(T>  (3.2) 

where  G(*)  is  the  mixing  distribution.  Because  (R^,  ...,  Rt>  are  con¬ 
ditionally  Independent  and  identically  distributed,  we  have 

P(N_  >  t)  -  /?[P(R/to.  <  c  |T)]CdG(T) 

R  0  0  a  ^  (3.3) 

>  [P(R/aQ  <  ca)]  -  P(Nr  >  t) 

* 

as  in  the  proof  of  Theorem  2.1,  where  N  has  the  distribution  G(t;a  ). 

R  R 

2 

Not  only  do  the  run-length  distributions  of  R  and  S  charts  depend  on 
the  underlying  spherical  process,  but  the  geometric  bounds  of  Theorem  3.2 
are  similarly  dependent  through  and  ctg.  For  these  bounds  to  be  useful 
it  thus  would  appear  that  the  form  of  the  spherical  process  must  be  known 
in  each  particular  case.  Our  next  developments  show  that  run-length  distri¬ 
butions  are  ordered  stochastically  whenever  the  underlying  distributions  are 
ordered  in  their  degrees  of  peakedness.  It  follows  that  geometric  bounds 
for  a  reference  distribution,  such  as  the  Gaussian  law,  then  apply  uniformly 

for  all  spherical  laws  more  peaked  than  the  reference  distribution. 

N 

Let  p(*)  and  v(-)  be  measures  on  R  ;  u(’)  is  said  to  be  more  peaked 

N  N 

about  9  e  R  than  v(*)  if,  for  every  compact  convex  set  E  c  R  symmetric 

about  9  under  reflection,  u(E)  >_  v(E)  (cf.  Sherman  (1955)).  If  u(*)  and 

v(*)  are  finite-dimensional  measures  characterizing  two  zero-mean  spherical 

processes,  then  the  u-process  is  said  to  be  more  peaked  about  zero  than 

the  v-process  if  U(*)  is  more  peaked  about  0  e  than  v(*)  for  every  N  <  °°. 

Our  basic  result  is  the  following. 

THEOREM  3.3.  Let  Fy(t;y),  with  paremeter  sequence  y  -  (y^,  y2>  •••),  be  the 
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2 

run-length  distribution  for  either  the  R  or  S  chart  of  a  drifting  process 
generated  from  a  spherical  process  having  the  measure  u(*).  Of  two  pro¬ 
cesses  having  the  measures  u(«)  and  v(*)>  suppose  u(*)  is  more  peaked  than 

2 

v(*).  Then  the  run-length  distributions  for  the  R  and  S  charts  from  the 
two  processes  are  ordered  stochastically  as 


F^Ctjy)  1  Fu(t;y) 


for  each  fixed  y  and  t  >  0. 


Proof.  Let  N  be  the  run  length  of  the  R  chart;  S  charts  are  treated 

. . .  K 

similarly.  Let  N  *  tn  and,  as  in  the  developments  preceding  Theorem  3.1, 

consider  {X^,  X^,  ....  Xt>  and  {Y^,  Y^,  •••,  Yfc>  such  that  Y^  »  y'^X^,  and 

let  R(Y.)  be  the  range  function  R(Y, )  •  Y. .  ,  -  Y.  ....  As  the  sample  range 
-l  -i  ivn;  1UJ 

and  sample  variance  are  invariant  under  translation,  we  may  take  the  process 
to  be  symmetric  about  zero.  Let 

At  -  (Y  e  RN|R(Yi)  <^1;  i  -  1,  2,  ...,  t 
and  A  -  A^  Under  the  typical  measure  ui  ( • )  on  R*1  and  the  induced  measure 
P  (• )  on  R*  we  infer  that 


P  (N  >  t)  -  u(A). 

U) 

Because  R(Y^) ,  the  sum  of  convex  functions,  i9  convex  and  symmetric,  it 
follows  that  A^  is  a  cylinder  set  in  R1*  symmetric  about  0  and  having  compact 
convex  sections.  Therefore  A  is  a  compact  convex  symmetric  set  and,  under 
the  hypothesis  that  u(')  is  more  peaked  about  0  e  /  than  v(*)»  it  follows 
that  u (A)  m(A),  i.e.. 


V*  ’  £)  *  C) 

which  is  equivalent  to  the  assertion  of  the  theorem. 


4.  Numerical  Comparisons.  We  study  numerically  some  bounds  developed 

in  earlier  sections,  first  considering  stochastic  bounds  on  F  (t)  and  Fc(t), 

R  b 
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as  given  in  Theorem  2.1,  when  Che  underlying  distributions  are  Gaussian. 

We  treat  the  case  that  n  •  10  observations  are  taken  on  each  occasion, 

2  2 
Oq  having  been  estimated  using  the  sample  variance,  SQm>  of  a  variable 

number,  m,  of  observations  in  a  base-line  period.  To  monitor  at  the 

level  ct  *  0.05,  we  choose  the  critical  values  c  (m)  and  d  (m)  from  standard 

a  a 

2  2 

tables  such  that  P(R. /S»  <  c  (m))  ■  0.95  ■  P(S. /Sn  <  d  (m)).  We  assume 

i  0m  —  a  i  Um  —  a 

2  2 

that  the  process  is  under  control,  i.e.,  ■  a q  for  all  i.  The  run-length 

2 

distributions  for  the  R  and  S  charts  are  found  by  numerical  integration 

from  normal-theory  versions  of  expressions  (2.1)  and  (2.2)  using  an  IMSL 

(1979)  subroutine.  Results  of  these  computations  are  tabulated  in  Table  1 

2 

for  R  charts  and  in  Table  2  for  S  charts  for  various  choices  of  t  and  m. 

We  report  results  up  to  t  «  100  and  m  ■  40,  as  computational  errors  are 

greater  for  large  values  of  t  and  m  owing  to  the  configuration  of  the  system 

2 

used.  The  entry  m  *  ®  in  Tables  1  and  2  corresponds  to  the  case  that  o^  is 

known,  the  base  period  remains  unobserved,  and  each  run-length  distribution 

is  geometric  having  the  parameter  a  *  0.05.  The  distribution  G(t;a)  also  is 

the  limit  of  F  (t)  and  F  (t)  as  m  -*■  «,  whatever  be  the  choice  for  a. 

R  S 

These  tabulations  suggest  a  further  ordering,  namely,  that  N  and  N 

R  b 

decrease  stochastically  as  m  increases,  reflecting  successively  weaker 

2  2  2  2 

dependencies  among  the  ratios  {R./Srt  ,  ...,  R  /S.  }  and  { S , / S_  ,  ....  S  /Sn  }. 

l  um  t  um  i  um  t  Dm 

The  geometric  approximation  is  seen  to  be  quite  close  even  for  m  *  40.  In 

addition.  Tables  1  and  2  suggest  that  the  R  chart  may  signal  more  frequently 
2 

than  the  S  chart  when  the  process  is  in  control,  its  variance  is  estimated 
2 

using  SQm  from  a  base  period,  and  both  charts  operate  at  the  level  a  *  0.05. 

We  next  study  examples  of  the  drifting  processes  considered  in  Section  2, 
specifically,  the  Gaussian  case  in  which  successive  samples  are  Independent. 
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We  suppose  chat  n  ■  10  observations  are  drawn  on  each  occasion,  including 

2  2 

the  base  period  which  yields  the  estimate  S.  for  o-.  Again  c  and  d  are 

0  0  a  ct 

taken  from  standard  tables  with  a  m  0.05.  To  study  the  bounds  of  Theorem 

2.3  numerically,  we  consider  the  run-length  distributions  F  (t;y)  and  F  (t;y) 

2  2 

associated  with  parameter  sequences  of  the  special  type  y(£)  ■  (5  ,  £  ,  ...), 
where  £  ■  ci/°0"  tables  3  and  4  give  probabilities  from  these  distributions 
for  £  taking  values  in  the  set  {1.0,  1.2,  1.5,  2.0},  entries  having  been 
computed  by  numerical  integration  using  an  IMSL  (1979)  subroutine.  These 
values  in  part  are  presented  graphically  in  Figures  1  and  2,  which  provide 
envelopes  for  the  run-length  distributions  of  certain  families  of  drifting 
processes*  For  example,  all  drifting  processes  which  satisfy  o^/o^  ^1.2 
for  all  i  will  have  run-length  distributions  lying  in  the  region  bounded  by 
the  curves  labeled  £  »  1.0  and  £  -  1.2.  Note  that,  while  the  run-length 
distributions  for  drifting  processes  depend  on  the  path  of  the  drift,  namely, 

Y  *  (y^,  Y2»  •••),  the  envelopes  depend  only  on  the  maximum  and  minimum  drifts 
of  the  process. 

Turning  to  the  spherical  processes  discussed  in  Section  3,  we  note  by 

Theorem  3.2  that  geometric  bounds  are  available  for  the  run- length  distribu- 

2 

tions  of  the  R  and  S  charts,  the  parameters  and  a  depending  on  the 

K  o 

particular  process.  It  is  instructive  to  study  disturbances  in  a  and  a 

R  S 

owing  to  misspecification  of  the  underlying  process.  Here  we  compare  the 
actual  values  with  their  nominal  values,  the  control  limits  ca  and  da  having 
been  set  under  the  mistaken  assumption  that  the  process  is  Gaussian,  when, 
in  fact,  successive  observations  are  generated  from  a  spherical  Cauchy  pro¬ 
cess.  These  comparisons  are  given  in  Table  5  for  R  charts  and  in  Table  6  for 
2 

S  charts  using  samples  of  sizes  n  ■  3,  5,  and  10.  Table  5  was  generated  by 


17 


IJ 


simulating  P(R/Oq  <_ c^) ,  c a  being  the  normal-theory  control  limit,  using 

100,000  replications  for  n  ■  3  and  n  ■  5  and  using  10,000  replications 

2  2 

for  n  ■  10.  Table  6  was  constructed  on  computing  P(S  /o^  <_  d^) ,  being 

2  2 

the  normal-theory  control  limit,  using  the  fact  that  S  / a ^  has  an  F-dis- 
tribution  with  n  -  1  and  1  degrees  of  freedom  in  the  case  of  a  spherical 
Cauchy  process.  The  evidence  at  hand  indicates  that  disturbances  in 
and  a,,  become  larger  as  the  sample  size  increases  under  the  particular 
type  of  misspecification  considered  here. 
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TABLE  1.  Run-length  distributions  of  R  charts  using  a  ®  0.05,  samples 


of  size  n  ■  10,  and  the  sample  variance  from  a  base  period  of  m  observa¬ 
tions,  under  normality  when  the  process  is  in  control. 

m 

P(Nr  >  5) 

P(Nr  >  10) 

P(Nr  >  20) 

P(Nr  >  50) 

P(N  >  100) 

K 

5 

0.90106 

0.85295 

0.79521 

0.74538 

0.73882 

10 

0.85438 

0.79618 

0.72168 

0.61366 

0.55701 

20 

0.82050 

0.72422 

0.60897 

0.45363 

0.33984 

40 

0.79660 

0.66909 

0.51149 

0.30223 

0.17738 

CO 

0.77378 

0.59873 

0.35848 

0.07694 

0.00592 

20 


2 

TABLE  2.  Run-length  distributions  of  S  charts  using  a  3  0.05,  samples 
of  size  n  *  10,  and  the  sample  variance  from  a  base  period  of  m  observa¬ 
tions,  under  normality  when  the  process  is  in  control. 


m 

P(Ng  >  5) 

P(Ng  >  10) 

P(Ng  >  20) 

P(NS  >  50) 

P(NS  » 

5 

0.90076 

0.87008 

0.83509 

0.80584 

0.79672 

10 

0.86113 

0.80833 

0.74950 

0.67011 

0.61025 

20 

0.82715 

0.74039 

0.63967 

0.50024 

0.40076 

40 

0.80168 

0.68204 

0.53659 

0.34162 

0.21937 

CO 

0.77378 

0.59873 

0.35848 

0.07694 

0.00592 

If 
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TABLE  3.  Run-length  distributions  for  the  R  charts  of  stationary  Gaussian 
processes  having  the  parameters  y(£)  =  (£2,  £2,  ...),  using  the  sample 
variance  in  a  base  period  of  m  =“10  observations  to  estimate  with 
a  =*  0.05  and  samples  of  size  n  *  10. 


£ 

P(Nr  >  5) 

P(Nr  >  10) 

P(Nr  >  20) 

P(Nr  >  50) 

o  1 
: 

W  1 

1.0 

0.85438 

0.79618 

0.72168 

0.61366 

0.55701 

1.2 

0.67862 

0.57321 

0.47783 

0.35281 

0.26283 

1.5 

0.38038 

0.25958 

0.16908 

0.09150 

0.05167 

2.0  0.08241 


0.03287 


0.01212 


0.00282 


0.00099 


TABLE  4.  Run-length  distributions  for  the  S  charts  of  stationary  Gaussian 
processes  having  the  parameters  Y(C)  *  ...)»  using  the  sample 

variance  in  a  base  period  of  m  =' 10  observations  to  estimate  o^, 
a  =*  0.05  and  samples  of  size  n  =  10. 


€ 

P(NS  >  5) 

P(Ng  >  10) 

P(Ng  >  20) 

P(Ng  >  50) 

P(NS  > 

1.0 

0.86113 

0.80833 

0.74950 

0.67011 

0.61025 

1.2 

0.68704 

0.59400 

0.50344 

0.39644 

0.32490 

1.5 

0.37996 

0.26980 

0.18574 

0.10983 

0.07304 

2.0  0.07654 


0.03249 


0.01315 


0.00386 


0.00151 


I 
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2  2 

TABLE  6.  Values  of  P(S  /Oq  <_  da)  for  spherical  Gaussian  and  Cauchy 
processes  using  samples  of  size  n. 


Gaussian 

Process 

n  =  3 

Cauchy  process 
n  *  5 

n  -  10 

0.900 

0.5778 

0.5130 

0.4538 

0.950 

0.6218 

0.5515 

0.4843 

0.975 

0.6546 

0.5814 

0.5089 

0.990 

0.6870 

0.6123 

0.5353 
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Figure  2.  Envelopes  for  Che  run  lengths  of  S  charts  of  drifting  Gaussian 
processes  in  terns  of  stationary  processes  having  the  parameters  7(5)  - 
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